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ABSTRACT 

We develop and apply a novel method of analysis to study the X-ray spectral 
properties of 908 QSOs in the Large Bright Quasar Survey (LBQS) that were observed 
during the soft X-ray ROSAT All-Sky Survey (RASS). Due to the relatively short 
(< 600 s) RASS exposure times, only 10% of the QSOs are detected in X-rays, so X- 
ray spectral model fits for individual QSOs are precluded by poor photon statistics. 
Spectral stacking provides effectively much more sensitive X-ray observations for an 
average QSO in bins of redshift, and for several classes of QSOs. We model the stacked 
X-ray spectra in a way that obviates both the distinction between detections and non- 
detections, and the need to stack together only those objects observed through similar 
Galactic column densities Nh- In appHcation to the LBQS/RASS sample, we marginally 
confirm a flattening of the X-ray spectral index of QSOs toward higher redshifts. Radio- 
loud QSOs show flatter X-ray photon indices, in agreement with previous work. We show 
that the apparent flattening of the photon index with redshift is not due to an increasing 
fractional contribution from the radio-loud class of QSOs. 

Key words: galaxies: active - quasars: general - X-rays: galaxies 
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1 INTRODUCTION 

Although quasi-stellar objects (hereafter, QSOs) radiate similar amounts of energy in the soft X-ray and optical bandpasses, the 
soft X-ray properties of these objects are by comparison poorly understood. Since the vast majority of soft X-ray observations 
of QSOs are broadband measurements of relatively low spectral resolution, the overall shape of the continuum is critical to 
calculation of X-ray flux and luminosity through counts-to-flux conversion factors and the cosmological ^-correction. These 
data are necessary for calculation of the intrinsic emission properties of QSOs, their X-ray luminosity function and its evolution 
with look-back time, and the contribution of QSOs to the diffuse X-ray background (XRB). The X-ray spectral slopes, and 
their connection to spectral slopes in other wavebands help to constrain the diverse physical models proposed to explain the 
continuum and line emission from the active galactic nucleus. 

The comparatively low energy resolution of most X-ray satellites has limited the analysis of QSO X-ray spectra to simple 
models (usually power law plus absorption). Low photon flux has further limited these studies mostly to apparently bright, 
low-redshift QSOs. As a result, the QSOs with well-studied X-ray spectra have necessarily formed a highly incomplete sample, 
biased by a variety of poorly-understood selection effects. Previous efforts to minimize selection effects have concentrated on 
complete samples of low-redshift QSOs observed in targeted exposures or all-sky surveys (e.g., Laor et al. 1994, Walter k. 
Fink 1993, respectively), or on deep pointings toward high redshift QSOs (e.g., Wilkes et al. 1992, Elvis et al. 1994a). 

To ease the deconvolution of selection effects from intrinsic properties, and to enable robust statistical analysis, a large 
sample of QSOs selected using uniform, quantiflable criteria and with detailed multiwavelength information available is 
desirable. The Large Bright Quasar Survey (Hewett, Foltz & Chaffee 1994, Morris et al. 1991, Chaffee et al. 1991, Hewett et 
al. 1991, Foltz et al. 1989, Foltz et al. 1987) is a large, uniformly-selected sample of QSOs covering a wide range of redshifts. 
LBQS QSO candidates were selected from photographic and objective prism plates, using a combination of quantiflable 
selection techniques sensitive to colours, to strong emission lines, redshifted absorption features or continuum breaks. As 
described in the series of LBQS papers, the technique appears to be highly efficient at flnding QSOs with 0.2 < z < 3.3, a 



4 N. Schartel et al. 

significantly broader range than past work. The entire LBQS optical sample consists of 1056 QSOs with 16.0 < Bj < 18.9 
over more than 450 deg^ of the sky. Measurement of the ensemble X-ray properties of QSOs to redshifts greater than a few 
tenths has so far been thwarted by the cost and difficulty of making numerous deep targeted X-ray exposures. Two approaches 
are available to circumvent these problems using a large, homogeneous sample such as the LBQS. 

One approach to determine the soft X-ray properties of such a sample is to search for deep X-ray images that serendip- 
itously include sample QSOs. A random subsample of 146 LBQS QSOs included in Einstein IPC pointings was studied by 
Margon et al. (1992). The small number of counts for most of the QSOs in the LBQS/IPC overlap sample, along with uncer- 
tainties in the response corrections for off-axis sources, preclude the derivation of detailed spectral constraints for individual 
QSOs. They co-added X-ray images of all QSOs in similar ranges of Galactic hydrogen column density N§"' maintaining all 
the spectral and background information. The properties of the 'stacked' QSO are then analysed as if it were a single object, 
but measurements instead provide average values characterizing the subsample of co-added QSOs. 

In a second approach, we study a much larger fraction of the LBQS sample now available for study using the first all-sky 
soft X-ray survey. The Position Sensitive Proportional Counter (PSPC, Pfefferman et al. 1987) onboard the ROSAT X-ray 
satellite (Triimper 1983, Aschenbach 1988) provides soft X-ray information for several thousand previously cataloged QSOs 
in the ROSAT All-Sky Survey (hereafter, RASS). The X-ray and optical properties of an LBQS/RASS overlap sample of 908 
QSOs have been described and analysed by Green et al. (1995). Due to the short (< 600 s) RASS exposure times, only 10% of 
the QSOs are individually detected in X-rays. However, by stacking X-ray counts, they obtain efi^ectively much more sensitive 
observations for an average QSO in bins of redshift or luminosity, and for several classes of QSOs. In the current work, we 
use stacking to study the soft X-ray spectral properties of several subsamples of the LBQS/RASS QSOs. 

Since the RASS exposure times for LBQS QSOs are short, even stacking in bins of N§"' often does not yield a sufficient 
number of counts to allow X-ray spectral analysis. Consequently, we have developed and introduce here a novel technique 
that permits us to co-add QSO X-ray spectral data, whether detected or not, and even when observed through substantially 
difi^erent Galactic column densities. In summary, we first stack together X-ray counts preserving the pulse height channel 
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(energy) information of each photon. The fit is performed using a multiplicative model that incorporates the absorption due 
to the Galactic column for each QSO, weighted according to its exposure time. This exposure-weighted absorption (hereafter, 
EWA) model is used throughout our analysis. To find the best-fit slope for a given spectral model, we fold the model through 
the PSPC response function for comparison to the data. If, after stacking, the signal-to- noise (S/N) ratio is still too low for a 
full least-squares fit to the observed X-ray spectrum, we use an adaptive hardness ratio (AHR) technique (see appendix § A) 
to derive the corresponding X-ray photon index. Our method is detailed in the following section (§ 2). We test the method 
with simulated high S/N data in § 3, and with low S/N data from the LBQS/RASS in § 4.1 Finally, we apply the EWA 
method to study the X-ray spectral properties of the LBQS QSOs in § 5 Readers most interested in the scientific results are 
directed to that section and the discussion following (§ 6). 

2 X-RAY SPECTRAL ANALYSIS 
2.1 The Spectral Stacking Method 

To generate the stacked spectrum from a sample of N sources, for simplicity we first describe the spectrum of a single source. 
As detailed in Green (1995), we used a source aperture radius of 180", optimized for optimal signal-to-noise of faint point 
sources in the RASS. A background annulus from 600 to 1200" is taken at each optical QSO position, from which unrelated 
detected sources are excised. In an exposure time t^, the number of counts received from the nth source in the ith pulse 
height channel is given by: 



where these variables represent for the source [background] of the nth source, the aperture area or observed solid angle 
[A^], the total number of photons Kn [Ln\, the vignetting correction factors V^^j. [Viti] for the A:th [Zth] photon, and the 
exposure time [f^]. The indices q, b and s pertain to the QSO, the background aperture, and the source aperture (total 
QSO plus background counts), respectively. PSPC background counts due to particles are not subject to vignetting as are 
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X-rays, but are correctly subtracted here for point sources. This is especially justified in the case of stacked RASS data, 
because in the RASS scanning mode, the source and the background are not observed at different off-axis angles, resulting in 
the same mean vignetting correction. 

The uncertainty in the net number of source counts Cf^ is 

\k = l " " 1=1 

The total exposure time is given by 

Af 
n = l 

where N is the number of sources. By summing over all N sources the counts in pulse height channels 8 - 240 (corresponding 
approximately to energies 0.08 to 2.4 keV) and dividing by the total exposure time T^, we get the stacked source spectrum, 

= {Ri, iig, ...ii24o}. The mean count rates for the individual channels for all N sources in the sample is thus 

N / K„ AS B ^-^^ \ 

R" - — ( X^V" - ^" v'' 1 

with error 

(N / K„ ^" 4s s 2 

E E E(tl) ("'■■' 
n=l \*:=1 (=1 

As an example, the stacked spectrum of all 72 detected QSOs in the LBQS/RASS sample with 0.2 < z < 1.4 is shown in 
Figure 1. 

2.2 Exposure- Weighted Absorption (EWA) Models 

Because the absorption of soft X-rays by material along the line of sight to a QSO is energy-dependent and a non-linear 
function of the intervening column density, mean spectral parameters from stacking have been determined either by stacking 
together only QSOs with similar N^"' (e.g., Margon et al. 1992) using 21cm measurements (e.g. Stark et al. 1992), or by 
stacking contours from fits to individual QSOs (Canizares &; White 1989). Here we outline a new method to derive the 
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mean spectral parameters for an ensemble of sources where the individual sources need not be detected, and may be observed 
through a range of absorbing columns. Although the method we outline is suitable for a variety of X-ray spectral models (e.g., 
Raymond-Smith, thermal bremsstrahlung, blackbody) we assume here a power law index, which generally provides the best 
fit to the X-ray spectra of QSOs. With this assumption, the model spectrum of the n-th source is given by: 



where /„ is the photon flux density of the n-th model spectrum in photons cm^^ s keV^^, Eq is the reference energy in 
keV, In is the model normalization (the photon flux density at Eq in photons cm~^ s ~^ keV~^), (F) is the mean photon 
index, {NH)n is the equivalent hydrogen column density of the n-th model spectrum (in cm~^), and a (E) is the photoelectric 
cross-section following Morrison &; McCammon (1983). We stress that the method is only rigorously correct if each input 
spectrum has a photon index equal to the mean (F). However, within the uncertainties, a mean index representative of the 
sample mean will result 

The model for the stacked spectrum is given by the sum of all N model spectra weighted by their exposure time: 



where F is the photon flux density of the model for the stacked spectrum in photons cm ^ s ^ keV ^ . Using the deflnition 
of the mean of the normalization (7), 




(1) 




dE 



(2) 




one can write for I, 



In = {I) + 51, 



n 



(3) 



where 



5In = 



(4) 



n=l 



Substituting equation 3 into equation 2 gives 
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= W(t Ite-'"-'-'"') dE 

If Sin is independent of {NH)n and i^, then, using equation 4, 

1^1^ <5/„ e-(^«)-(^) j «0 (5) 
one gets the model of the stacked spectrum: 

F{E,{NH)n,In,{r),Eo,N) dE^ ^^-(NhU<^(e)^ . (6) 

The model of the stacked spectrum is therefore given by the convolution of the exposure-weighted absorption with a 
power law whose normalization is the mean of the individual source normalizations. 

We note that to meet the requirement for EWA models that In be independent of N§"' and tn (see equation 5), samples 
which are X-ray-selected must be treated with caution. For example, a sample consisting of sources required to have more than 
some number of X-ray counts will show a correlation between /„ and N§°', since only those with larger intrinsic luminosities 
will meet the selection criterion. Such a problem may be partially alleviated by the randomizing factor of distance, but is best 
avoided by using samples selected in other wavebands, such as the LBQS. 



3 TESTS WITH HIGH SIGNAL TO NOISE SIMULATED SPECTRA 
3.1 Tested Models 

Before turning to the relatively low S/N RASS observations of the LBQS QSOs, we first test our technique using high S/N 
simulated data. We simulate the ROSAT PSPC spectra of 4 individual QSOs assuming a variety of power law spectral slopes 
and intervening absorbing columns. The simulated spectra are then stacked as described above. Finally, we fit the stacked 
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spectrum with 4 power law models to see how well each reproduces the original input spectral parameters. First, we test a 

power law with 3 free parameters (normalization J, Nh and F, the 'free TVh' fit). Second, we use a model with absorbing column 

JV ^ 

fixed at the exposure-weighted mean Nh- By this, we mean equation 1, but with Nh replaced by (Nh)^^ = ^{NH)n- 

n=l 

Third, we test a model with the absorbing column fixed at the counts- weighted mean Nh- By this, we mean equation 1, but 
with Nh replaced by {Nh)'^^ = Yj %I-{NhU. The latter model can be used only when the true source counts are known 

n=l 

for each QSO, i.e., when every object in the sample is detected. Finally, we test our own exposure- weighted absorption (EWA) 
model, described above, which can be used irrespective of sample detection fraction. 

3.2 Simulated Samples 

We are most often interested in investigating the X-ray spectra of samples with similar non-X-ray parameters. For example, we 
investigate below the dependence of the average QSO X-ray power law slope on radio-loudness and redshift. The investigation 
of spectral evolution with look-back time requires binning in redshift. (Such binning also prevents substantial degradation in 
spectral resolution from co-adding a wide range of restframe energies.) Since in any magnitude-limited sample, redshift and 
flux are strongly correlated, we infer that QSOs in a redshift bin will tend to have similar X-ray fluxes. Since it is generally 
known that lopt and Ix are strongly correlated, this clustering in X-ray flux will obtain also for optically-selected sample 
such as the LBQS (Green et al. 1995). For the particular case of the RASS, the exposure times are also similar. For these 
reasons, we have chosen our primary simulated samples to reflect this by requiring fixed exposure time and intrinsic flux. 

To test how well stacking and EWA models reproduce the simulated parameters in a variety of situations, we generate 
3 representative samples of 4 QSOs, each with fixed intrinsic fiux and exposure time t. Sample A [t = 500 s) shares a single 
power law slope (F = —2.5) but covers a range of Nh- Sample B QSOs {t = 750 s) have Nh = 2.5 x 10^° cm~^ and a range 
of F. Sample C (t = 750 s) covers a range of both F and Nh- We then fit the individual spectra. Due to the intentionally high 
S/N, the individual fits reproduce the model parameters to within < 1%. After fitting, the number of counts of each spectrum 
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Table 1. Input parameters for simulated QSO spectra. 



Sample A Sample B 

(r = -2.5, vary TVh) {Nh ^2.5'' , vary T) 
iV° Counts r Counts 



Ma 



Sample C 
(vary Nh , vary T) 

r Counts 



1.0 378693 

3.0 174568 

5.0 119193 

7.0 96498 



-2.00 218798 

-2.25 252554 

-2.75 370993 

-3.00 472399 



1.0 
7.0 
5.0 
3.0 



-2.00 
-2.25 
-2.75 
-3.00 



344151 
138599 
199497 
390066 



" Intervening column density in units lO^^cm ^. 

is adjusted so that each source has an intrinsic monochromatic flux 7 at 1 keV of 1.0 photon cm s^i keY-\ In Table 1, we 
show the model parameters and resulting counts for the 12 QSOs in the 3 samples. 

After stacking the spectra in each sample according to the prescription in § 2.1, we then flt each stacked spectrum using 
4 power law models (with free Nh, with exposure- and counts- weighted {Nh), and with EWA). Figure 2. illustrates these flts, 
with the best fit parameters for each sample/model combination listed in Table 2. 

For sample A, EWA clearly provides the best fit method. EWA models are thus particularly well-suited to modeling stacked 
spectra of sources expected to have similar intrinsic spectral properties, when observed at many diflFerent sky positions or 
through patchy absorption. 

For sample B, EWA and {Nh) models yield identical results as expected, since the effective absorption in both models 
is identical. The model with Nh free inaccurately reproduces the input Nh, but appears to reproduce the input F, and I 
parameters somewhat more accurately. Sample B might be approximated by a sample of QSOs observed in a single deep 
X-ray image. In this case, we recommend a model with Nh free to best measure F, but only if there are adequate counts for 
a two-parameter fit. 
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Table 2. Power law fits to stacked simulated samples. 



Sample 


Absorption 




Best Fit Parameter" 






Model 


Nh 
(1020 cm-2) 


/ 

( cm-2 s-i keV-i) 


r 


xVd.o.f. 


A 


free Nh 


1.84 ± 0.02 


n Q'ifi -t- n nn'i 


— z.zoo ziz u.uuy 


55.5 / 25 


A 




4.0 


n Qf;7 -t- n nn'i 


Q n9fi -t- n nn'i 

— o.uzu zn u.uuo 


800.2 / 26 


A 




2.83 


n Qf;7 -t- n nn'i 

u.yoi zn u.uuo 


— z.uio zn u.uuo 


242.0 / 26 


A 


EWA 




1 nni -L n nnQ 
j-.uuj. zn u.uuo 


_9 c:n9 _|_ n 009 
z.ouz zn u.uuz 


18.5 / 26 


B 


free Nh 


2.29 ± 0.02 


1.004 zt 0.003 


-2.521 zt 0.007 


7.1 / 25 


B 




2.5 


1.008 zt 0.002 


-2.600 zt 0.002 


22.0 / 26 


B 




2.5 


1.008 zt 0.002 


-2.600 zt 0.002 


22.0 / 26 


B 


EWA 




1.008 zt 0.002 


-2.600 zt 0.002 


22.0 / 26 


C 


free Njj 


2.43 ± 0.02 


0.965 zt 0.003 


-2.364 zt 0.007 


12.8 / 25 


C 




4.0 


0.985 zt 0.002 


-2.907 zt 0.002 


545.0 / 26 


C 




3.25 


0.981 zt 0.002 


-2.651 zt 0.002 


172.1 / 26 


C 


EWA 




1.013 zt 0.002 


-2.387 zt 0.002 


112.8 / 26 



" Errors are taken from the covariance matrix of the fit and correspond to the 68.3 % confidence level. 

In the most general case, Sample C, where both objects and sight-lines may have a dispersion in properties, EWA does 
not provide the best fit as judged by values, but most accurately retrieves the true (F). If all sources are detected, the 
intrinsic properties of the sample may be adequately understood by simply assuming a count-weighted mean Nh. But the 
^■p^cw luQi^el cannot be used when some sources have only X-ray fiux upper limits, as for the LBQS/RASS sample. 

We have also tested an analogous set of samples, where individual QSO spectra may have different exposure times, but 
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are generated to have the same total counts and flux normalization. For such samples, EWA again provides best-flt spectral 
parameters closest to the input simulated spectra. 

Overall, the EWA model performs as well or better than the other 3 power-law models tested for these samples, and is 
clearly superior when considering sources with few counts. We expect that strong average spectral constraints for a variety of 
object classes (not just QSOs) might be achieved by stacking spectra of samples observed in the large number of Einstein or 
ROSAT deep pointed observations now becoming available in data archives, and by using EWA models to correctly account 
for intervening absorption. 

3.3 Absorption: Intrinsic or Galactic? 

For no QSO in the LBQS/RASS sample are there sufficient counts for adequate power law spectral fits to constrain both 
r and the total line-of-sight column, N*jJ*=N'j?*''+ Ng"'. Thus, we (and many other workers) proceed with the assumption 
that the intrinsic column is negligible in the aggregate. The majority of optically-selected QSOs observed in soft X-rays to 
date indeed show little evidence for strong absorbing columns N^J*'' intrinsic to the QSO (e.g., Fiore et al. 1994, Walter &; 
Fink 1993, Canizares & White 1989, Wilkes & Elvis 1987). We verify this for the LBQS/RASS sample in § 4.1, and hereafter 
perform all model fits with absorption fixed toward individual sources at the Galactic value (taken from Stark et al. 1992), 
leaving two free parameters, the photon index F, and the flux normalization I. 

4 TESTS WITH LOW SIGNAL TO NOISE DATA 
4.1 Detection status and absorption 

Green et al. (1995) describe criteria for detection of LBQS QSOs (or other faint point sources) observed in the RASS, 
optimized for maximum S/N. When stacking, we co-add X-ray counts from all QSOs regardless of their detection status, 
so we should derive similar spectral indices when stacking together samples of either detected or non-detected QSOs, if 
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the intrinsic spectral behavior does not depend strongly on the measured count rate. Furthermore, if the derived spectral 
parameters differ for subsamples of QSOs observed through different N§"', this would suggest that fitting with an EWA model 
does not correctly treat the effects of absorption. 

As a check on our technique, we thus divide the LBQS/RASS sample into subsamples, one with 1.5 < N§'^' < 2.5 (in 
units of 10^° cm~^), a second with 2.5 < N^"' < 3.5, and a third containing QSOs regardless of N§"'. We further divide these 
according to detection status before stacking. Table 3 shows for each of these subsamples the number of QSOs, the range and 
mean of the Galactic column density N^"', the stacked count rate (counts per 1000 s), and the mean photon index from the 
best-fit EWA model (see § 2.2). Since there may be an intrinsic evolution of the mean X-ray photon index (F) with redshift, 
and since higher redshift QSOs are much less likely to be detected, we include only QSOs with 0.2 < z < 1.0. Figure 1. shows 
the spectrum that results from stacking all 72 detected QSOs in this redshift range. Because the total number of photons in 
the stacked spectra of non-detected QSOs is small, we determine the spectral parameters of these using the adaptive hardness 
ratio technique (see the appendix). 

We find that the best-fit photon index is, within the uncertainties, independent of either the division into absorption 
groups or the classification by detection status. A comparison of photon indices and their errors in Table 3 shows that the 
EWA method yields consistent photon indices across bins of Galactic column density (e.g., for the classes 1, 2, 4, and 5). A 
comparison between the mean photon indices of the detections, 2.65 ±0.07 (class 3), and the non-detections, 2.56 ±0.11 (class 
6), also shows no significant difference. As expected, our results are thus independent of count rate and intervening Galactic 
absorption. The EWA model can be rigorously applied to the stacked spectra of QSOs observed with a wide range of exposure 
times or through a wide range of Galactic column. 

We also performed fits to the stacked spectra including the possibility of absorption intrinsic to the QSOs. To do this, we 
convolve a power law with intrinsic absorption with the exposure-weighted Galactic absorption. The relatively small number 
of counts yields poor statistics, but fits were consistent at the la level for all of the above classes having no intrinsic equivalent 
hydrogen column density (i.e., N^* fa 7V^"'), justifying our explicit assumption of this in following analysis. 
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Table 3. Detection status and Galactic absorption. 



Class 
No. 


No. of QSOs" 


]\jGal b 

Range Mean 


Count rate 
(per ksec) 


(r>- 

(EWA) 


Xred. 1 ^■°-/- 






Detections 








1 


37 


1.5 - 2.5 1.99 


73.5 ±2.1 


2.68±0.10 


1.74 / 12 


2 


26 


2.5 - 3.5 3.03 


59.3 ±3.2 


2.61±0.13 


1.23 / 11 


3 


72 


all 2.53 


66.9 ±2.0 


2.65±0.07 


1.06 / 15 






Non-detections'* 








4 


129 


1.5 - 2.5 2.05 


10.5 ± 1.1 


2.62±0.16 


AHR 


5 


89 


2.5 - 3.5 3.05 


9.4 ± 1.2 


2.46±0.20 


AHR 


6 


275 


all 2.81 


9.6 ± 0.7 


2.56±0.11 


AHR 



" Only QSOs with z < I included in this analysis. 
* Galactic column density in units lO^^cm^^. 

Errors are taken from the covariance matrix of the fit and correspond to the 68.3 % confidence level. 

Photon indices for non-detections calculated using the adaptive hardness ratio (AHR) technique. 

5 SPECTRAL CHARACTERISTICS OF THE LBQS/RASS SAMPLE 
5.1 Mean photon indices for the LBQS: RASS vs. Einstein 

Stacking the spectra of the entire LBQS/RASS sample of 908 QSOs achieves a single spectrum with more than 6000 counts. 
Were so many counts available for a single QSO, or for an ensemble of QSOs in a very narrow redshift range, somewhat more 
detailed spectral analysis would be possible. The stacked spectrum of a sample such as the LBQS encompasses a wide range 
of redshifts, blurring out many spectral details. Since we are investigating assumed power law slopes that by their nature are 
conserved with respect to redshift, our derived mean X-ray photon index (r) should nevertheless adequately represent that 
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of the mean QSO in the LBQS. We thus perform a standard fit to the stacked spectrum, assuming no intrinsic absorption. 
The result (with its error) is (F) = 2.53 ± 0.07 (listed in the last row of Table 4). This is consistent with the mean (with its 
dispersion) of (T) = 2.5 ± 0.5 for a sample of 58 low redshift ((z) = 0.08) AGN measured in the RASS by Walter & Fink 
(1993). An initial PSPC spectral survey of 10 PG quasars by Laor et al. (1994) yields (F) = 2.5 ± 0.4, also consistent with 
our much larger sample. 

We may also compare our derived (F) for the entire LBQS/RASS with that of the serendipitous Einstein sample (the 
LBQS/IPC sample; Margon et al. 1992), since the latter includes a random subsample of the former. This is confirmed by 
the similarity of the redshift distributions for the LBQS/RASS sample (see Table 4) to that of the Einstein sample, where 
0.2 < z < 3.1, and {z) = 1.33. The mean photon index of the LBQS/RASS sample, Trass = 2.53 ±0.07, is measured between 
0.08 keV and 2.4 keV, while the LBQS/IPC analysis yields F/pc = 2.3 ± 0.2, measured between 0.16 keV - 3.5 keV (both 
bandpasses in the observer's frame). We thus find no significant difference between the mean photon index of the LBQS QSOs 
using Einstein and RASS spectra, but the combined uncertainty is large. 

A variety of papers (e.g., Schartel, Walter, & Fink 1993, Brunner et al. 1992, Laor et al. 1994) have found steeper (softer) 
power law slopes from ROSAT PSPC than Einstein IPC data. Fiore et al. (1994) demonstrate that fits to IPC spectra of 
QSOs yield indices lower on average by 0.5 — 0.7 compared to PSPC spectra even for the same objects. They argue that some 
of this discrepancy may be due to calibration errors in the response matrices of one or both of the IPC and PSPC, which may 
weaken the the significance of direct PSPC/IPC comparisons. 

5.2 Redshift subsamples 

We have derived the photon index for stacked subsamples of LBQS/RASS QSOs divided into 4 bins of redshift, with bin 
edges z = 0.2, 0.4, 0.6, 1.0, and 3.3. We chose these particular bins to optimize comparison with results for RQ and RL QSOs 
in § 5.3, but we note that our results are not significantly affected by the exact choice of redshift bin boundaries. The mean 
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Table 4. Redshift. 



Class 


No. of QSOs 


Redshift 


^total 


Count rate 




T GOi. ' 


No. 




Range 


Mean 


(ksec) 


(per ksec) 


(EWA) 




1 


85 


0.2 - 0.4 


0.32 


47.9 


34.4 ± 1.5 


2.51 ± 0.11 


1.16 / 10 


2 


89 


0.4 - 0.6 


0.50 


51.3 


23.3 ± 1.4 


2.74 ± 0.14 


0.68 / 8 


3 


174 


0.6 - 1.0 


0.80 


96.5 


15.1 ± 0.9 


2.78 ± 0.15 


0.88 / 4 


4 


560 


1.0 - 3.3 


1.72 


316.7 


5.8 ± 0.5 


2.28 ± 0.18 


0.98 / 4 


5 


348 


0.2 - 1.0 


0.61 


195.7 


22.0 ± 0.7 


2.66 ± 0.07 


AHR 


6 


908 


0.2 - 3.3 


1.29 


512.4 


12.0 ± 0.4 


2.53 ± 0.07 


AHR 



" Errors are taken from the covariance matrix of the fit and correspond to the 68.3 % confidence level. 

photon indices by redshift are shown in Table 4, along with the overall mean photon index from EWA models for the entire 
LBQS/RASS sample. 

Figure 3. shows the mean ROSAT photon index fits to stacked spectra from our four redshift subsamples (filled circles 
with solid bars). The width of the horizontal bars represents the redshift bin range for stacking. The mean photon index for 
the entire LBQS/RASS sample is also plotted (filled circle, dashed bars). The largest observed difference is a decrease at the 
« 2(T level between the photon indices for classes 3 and 4, with mean redshifts {z) = 0.8 and 1.9, respectively. To reduce the 
uncertainty in the mean low redshift index, we also tried a broader low-redshift subsample (0.2 < z < 1.0; class 5), for which 
we still find a difference of similar significance. 

The decrease in measured photon index with redshift is similar to that of Stewart et al. (1994) from PSPC spectra in 
5 deep ROSAT fields, containing about 100 UVX QSOs (Boyle et al. 1990). As Stewart et al. pointed out, the power-law 
component of QSOs to ziS 3 does not attain the medium-energy slope of the cosmic XRB over the 2-10 keV range (F = 1.4; 
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see Holt 1992 for a review), even though the rest-frame bandpass is similar. Thus QSOs cannot form 100% of the XRB. Auto- 
correlation of the XRB also suggests that if the remainder is to come substantially from point sources, they must have higher 
surface densities than QSOs (Georgantopolous 1993). Narrow emission-line galaxies may provide both the harder spectra and 
higher surface densities required (e.g.. Shanks et al. 1991, Jones et al. 1995). 

From a large sample of serendipitous detected sources in archived ROSAT PSPC pointings, Vikhlinin et al. (1995) 
present a study of X-ray spectral slope as a function of source flux. They flnd that, after some corrections for Galactic stellar 
contamination, the mean photon index decreases monotonically toward lower fluxes. A direct extrapolation of that trend to 
fluxes of about 10~^^erg cm~^ s~^ is consistent with the slope of the cosmic XRB, and includes all point-source contributors. 
Since Vikhlinin et al. (1995) invoke X-ray selection, however, we suspect that (the more X-ray luminous) RL QSOs may make 
a larger fractional contribution to their derived mean photon index (and also to the XRB), which in turn would yield lower 
values of (F). Even though the LBQS/RASS sample is optically-selected, we conflrm the hardening trend observed by Vikhlinin 
et al. in the flux decade where our samples overlap, between 10^^^ and 10^^^ erg cm^^ s^^. (The counts-to-flux conversion 
factor for our redshift subsamples in Table 4 is about 2.2 x 10~^^ cts~^ erg cm~^.) The large, deep sample of Vikhlinin et al. 
(1995) is complementary to the current study, and is clearly an important contribution to the ongoing debate as to the origin 
of the cosmic XRB. However, the importance of complete optical spectroscopic foUowup should not be underestimated. 

5.3 The radio-loud and the radio-quiet subsample 

Optically-selected samples of QSOs appear to be separated into two distinct populations, radio-loud (RL) and radio-quiet 
(RQ; e.g., Kellerman et al. 1989). Although these populations have quite similar optical properties, the X-ray emission from 
RL QSOs is now generally accepted to be stronger than that from RQ QSOs (e.g., Avni &; Tananbaum 1986, Worrall et 
al. 1987). In addition, there appears to be a signiflcant difference between the mean X-ray photon indices of heterogeneous 
samples of radio-loud and radio-quiet QSOs (e.g., Shastri et al. 1993, Molendi et al. 1992, Wilkes & Elvis 1987). We are now 
able to examine the ROSAT X-ray spectral properties of a substantial fraction of the LBQS sample, divided by radio-loudness. 
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Table 5. Radio loudness. 



Class 


No. of QSOs 




^total 


Count rate 


(r>- 










(ksec) 


(per ksec) 


(EWA) 




Radio-loud 


40 


1.22 


22.5 


21.5 ± 2.1 


2.15 ± 0.23 


0.41 / 3 


Radio-quiet 


147 


1.38 


85.4 


13.3 ± 1.0 


2.65 ± 0.18 


1.17 / 5 



" Errors are taken from the covariance matrix of the fit and correspond to the 68.3 % confidence level. 

Radio data are currently available for 226 of the QSOs in the LBQS/RASS sample from VLA observations at 8.4 GHz 
(Visnovsky et al. 1992, Hooper et al. 1995), and from 26 (much less sensitive) detections from the sky surveys of Condon et al. 
(1988, 1991) at 5 GHz, or Large et al. (1981) at 408 MHz. RASS data for this latter radio-loud sample was recently analysed by 
Brinkmann, Siebert, &; Boiler (1994) and Brinkmann &; Siebert (1995). We define an object as radio-loud if the log ratio of the 
emitted monochromatic luminosities at 8.4GHz and 4410 A, log R = log {Is. 4 ghz//^^^^^ ° ), is greater than unity, as described 
in Green et al. (1995). This definition yields 193 radio-quiet and 42 radio-loud QSOs. There is only 1 RL QSO below z = 0.4, 
compared with 20 RQ QSOs. Furthermore, z = 0.4 corresponds to the mean redshift of QSOs with Mb = —24: Hooper et al. 
(1995) and Delia Ceca et al. (1994) find that the fraction of RL QSOs decreases abruptly for absolute magnitudes fainter than 
this. There is also only 1 RL QSO above z = 2.4, compared to 26 RQ QSOs. For these reasons, to enable a fair comparison 
between the RL and RQ QSO samples, we restrict the redshift range of both samples to 0.4 < z < 2.4. This leaves 40 RL 
QSOs (with a mean and r.m.s. for the redshift distribution of (z) = 1.2 ± 0.5), and 147 RQ QSOs {{z) = 1.4 ± 0.6). 

Table 5 briefiy characterizes the RL and RQ QSO samples and power law model fits to their stacked spectra as follows: 
the mean redshift, the total exposure time, number of counts per 1000 s, and the photon index as derived assuming in the 
EWA models that Nj}* = N§'^' . The RQ QSOs show a steeper mean spectrum ((F) = 2.7±0.2) than the RL QSOs (2.2±0.2), 
but only at the 1.7a level. Although the RL QSOs in the LBQS/RASS show higher count rates (at 3.5(t, see Table 5), since RL 
QSOs are fewer in number, their total number of counts is smaller by a factor of about 2.5, resulting in larger uncertainties in 
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their mean photon index (F). Our results for RL and RQ QSO samples are consistent with the ROSAT PSPC photon indices 
found by other authors when similar redshift ranges are covered. These include studies from smaller and more heterogeneous 
samples of pointed observations (e.g., Fiore et al. 1994, Laor et al. 1994), as well as studies using RASS detections (e.g., 
Brunner et al. 1992, Schartel, Walter, & Fink 1993, Brinkmann, Siebert & Boiler 1994). Unlike the current work, all these 
studies suffer from X-ray selection effects, but they have the advantage of better overall photon statistics. 

5.4 QSOs with strong Fell emission 

Emission lines in QSOs are attributed to photoionization of gas by the high energy continuum near the QSO nucleus. The 
details of the process are not well-understood, but may be illuminated by comparison of emission line to continuum properties 
in samples of QSOs. Based on a heterogeneous sample of 18 objects with {z) = 0.15, Shastri et al. (1993) found that (radio- 
quiet) QSOs with strong Fell A4570 emission show softer (steeper) Einstein spectral slopes. We may contrast this result with 
X-ray spectral constraints on a subsample of the LBQS/RASS QSOs showing strong UV Fell emission. We use the subsample 
defined in Green et al. (1995), consisting of 23 QSOs in the LBQS/RASS sample flagged as having strong UV Fell emission. 
Although the Fell flag represents a subjective judgment, these QSOs are likely to be among those with the strongest iron 
emission in the LBQS sample. The redshift range most likely to lead to their Fell classification is 0.4 < z < 1.5, due to the 
iron feature under [Ne IV] A2423. We thus restrict the redshifts of both the Fell and a comparison non-Fell sample to within 
this range. The mean redshift of the Fell sample is {z) = 0.88. The comparison sample of 467 non-Fell QSOs has a similar 
mean redshift of 0.96. 

Using an EWA model and the AHR technique, we derive for Fell-strong QSOs a mean photon index (F) = 2.8 ± 0.2, 
consistent within the uncertainties with our result for the non-Fell sample, (F) = 2.7 ± 0.1. 

Our sample is probably more uniformly-selected that that of Shastri et al. (1993). Our sample has a wider range of 
redshifts, with a much larger mean redshift, and our classification is based on rest-frame UV, not optical, iron lines. The 
sample of Shastri et al. is limited to RQ QSOs only. Twelve of our 23 Fell-strong QSOs have radio data, of which 1148-0033 
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and 0107-0235 are radio-loud QSOs; only the latter is detected in the RASS. The exclusion of these 2 quasars from the stacked 
spectra does not significantly change our results. Neither would the exclusion of one or two more RL QSOs possibly included 
in the 11 QSOs without radio data, since no single QSO dominates the counts in the stacked spectrum. Since different energy 
ranges and different iron lines are sampled, we consider our null result to complement rather than contradict the results of 
Shastri et al. More detailed analysis, including correlations of Fell equivalent width and X-ray luminosity, are pending. 

6 DISCUSSION 

What might cause the apparent flattening (hardening) of the 0.08 keV - 2.4 keV X-ray slope we flnd toward higher redshifts? 
We consider briefly several possible explanations of the observed decrease in (r) that occurs somewhere between redshift 1 and 
2. We order this discussion by considering effects increasingly linked to the intrinsic physics of the QSO nuclear continuum. 

1) An increasing contribution from radio-loud QSOs. Since RL QSOs are known to have flatter X-ray spectra, and to be 
generally more luminous in X-rays than their RQ counterparts, the observed flattening of (F) toward higher redshift might 
simply be due to an increased contribution of X-ray counts from RL QSOs at higher redshift. Since stacking eliminates the 
distinction between detections and non-detections, this trend would not be caused by the larger detection fraction of the 
more X-ray-luminous RL QSOs. However, there is some evidence (e.g., Hooper et al. 1995) that RL QSOs may represent 
an increasing fraction of all QSOs at high redshift. Also, the study of Green et al. (1995) suggested for the LBQS/RASS 
sample that RQ QSOs show a relative decrease in X-ray luminosity with redshift not seen in RL QSOs. The fractional X-ray 
contribution of RL QSOs thus increases, possibly leading to the observed decrease of (F) with redshift. We test this possibility 
in two ways below. 

First, we examine whether trends in (F) exist within the RL and RQ populations separately. For a sample of RASS- 
detected RL QSOs, Schartel, Walter, & Fink (1993) indeed found a decrease in (F). In their highest redshift bin, {{z) « 1.7) 
they flnd (F) = 1.5 ± 0.2. It is intriguing that this value approaches that of the low-flux sources in Vikhlinin et al. (1995), and 
that of the XRB as well. Do RQ QSOs alone show a trend in (F) with redshift? Since the LBQS/RASS sample of RQ QSOs 
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is relatively large, we may stack their RASS spectra in 3 redshift bins, with edges at « = 0.2, 0.6, 1.0, 3.3. The photon indices 
for these 3 bins (with (z) = 0.40, 0.77, and 1.87) derived using the AHR technique are 2.7 ± 0.1, 3.2 ± 0.4, and 2.4 ± 0.3, 
respectively. Although there is some suggestion that the hardening occurs even for RQ QSOs alone, the apparent decrease in 
(r) is only at about the 1.4a level. 

Second, we estimate the fraction of X-ray photons in each redshift bin that are from RL QSOs. Let Rrl = NrlI{Nrl + 
Nrq), where Nrl {Nrq) is the number of photons attributable to RL (RQ) QSOs. Since the non-VLA subsamples contribute 
only detections, we must calculate Rrl for the VLA-observed subsample and assume that it applies to the overall LBQS/RASS. 
Where the largest change in (F) occurs, between QSO redshift classes of (z) = 0.8 and 1.7, Rrl decreases from 27.3% to 
20.9%. Such a decrease would be expected to produce an increase in (F), opposite to the trend observed. Thus, if this fraction 
for the VLA subsample is representative of the fraction for the whole LBQS/RASS sample in these redshift bins, the apparent 
flattening (hardening) of the X-ray spectra is not due to an increasing contribution of RL QSOs. 

2) Absorption over and above that due to our Galaxy. The small number of total counts necessitates our assumption 
here that iV^* « N^"' , leaving two fit parameters F, and /. This assumption is consistent with results from most previous 
X-ray studies of low redshift QSOs when fitting a single power law model. Exceptions occur in some sources with "warm 
absorbers" (e.g., Nandra &: Pounds 1992, Laor et al. 1994), and in a few intermediate to high-redshift RL QSOs (e.g., Elvis et 
al. 1994a,b). Consistent with these latter results, it is possible that the observed trend in (F) is due to absorption along the 
line of sight that increases toward higher redshift. Such absorption would reduce the observed soft X-ray emission, whereby 
our one-parameter fits here would yield a harder apparent photon index. If absorption is the sole cause of the observed trends, 
then it appears to become significant only for redshifts greater than about one. This would imply a strong evolution of the 
intrinsic absorption properties of QSOs with look-back time, or (due to the strong correlation between redshift and luminosity 
in the LBQS) with luminosity. Since strong X-ray absorption in high luminosity QSOs is unusual (e.g., Wilkes et al. 1987, 
Williams et al. 1992), intervening absorbers are more likely, as argued by Elvis et al. (1994a). The latter study was of RL 
QSOs only. Assuming that RL QSOs are probably giant ellipticals in clusters, those authors postulate that the absorption 
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is associated with cold intracluster gas (e.g., cooling flows). Since our overall counts are dominated by RQ QSOs, and since 
there is some weak additional evidence for a decrease in photon index with redshift for the RQ class alone, our work suggests 
that this apparent hardening is not solely due to associated absorption in the presumed cluster environments of RL QSOs. 
A hypothesis not inconsistent with our results is that of intervening absorbers such as damped Lya clouds at 0.3 <« < 1.9, 
(Elvis et al. 1994a, Vikhlinin et al. 1995). 

3) A soft excess due to a double power law. In many QSOs, at least two emission components are discernible in the X-rays: 
a flat nonthermal power law (F 1.7 — 2.0; e.g.. Turner &; Pounds 1989), and a steeper soft X-ray component (F 2.8; 
Fiore et al. 1993), with the spectral break between these two components somewhere between 0.3 and 1 keV (e.g., Jackson, 
Browne, &; Warwick 1993). A break energy in this range would produce a hardening of the RASS spectra in the redshift range 
1 < z < 2 as observed, since the soft component would shift out of the ROSAT band between those redshifts. 

4) The apparent hardening with redshift of the best-flt power law photon index (F) for QSOs in the LBQS/RASS sample 
might also be evidence for the evolution of the physical conditions in the QSO nucleus. For instance, the models of Haardt &; 
Maraschi (1993) predict that larger black hole masses yield softer X-ray spectra. This would be consistent with scenarios in 
which the black hole mass grows with time by accretion in long-lived QSOs, as discussed in Bechtold et al. (1994). 

In summary, we have studied here the soft X-ray spectral properties of QSOs in a large, homogeneous sample, the 
LBQS/RASS. We flnd that the QSOs in our sample with strong optical Fell emission have a mean X-ray spectral photon 
index (F) consistent with that of a control sample. We find some evidence for a difference in (F) between radio-loud and 
radio-quiet QSOs, and between QSOs at low and high redshift. Due to the moderate PSPC spectral resolution, and the 
relatively short exposure times of the RASS, the data here are not adequate to distinguish between the variety of proposed 
causes of the observed decrease in (F) with redshift. We have, however, demonstrated the utility of a novel spectral stacking 
technique which may improve this situation by facilitating rigorous analysis of co-added X-ray spectral data from samples of 
QSOs with deep pointed X-ray observations, even across wide regions of they sky. 
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APPENDIX A: FITTING THE MODEL TO THE DATA 

To determine the best-fit slope for a power law model, whether for an individual source spectrum or a stacked spectrum, 
we fold the model through the PSPC response function for comparison to the data. If, after stacking, the total counts are 
still too low for a full least-squares fit of a power law model to the observed X-ray spectrum, we can use the correspondence 
between hardness ratios and X-ray photon index to derive estimates for the latter quantity. We use two adaptive hardness 
ratios (AHRs): 
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HRl 



{B-A)/{A + B) 



and 



HR2 



{C-D)/{C + D) 



A, B, C, and D are the estimated number of (vignetting-corrected) source counts in the following channels: A : 8 ^ b, B : 
(6+1) — > 235, C : (6 + 1) — > c, and D : (c+l) —> 235. Here the energy of a given channel in keV is approximately (channel 
no.)/100. The channel numbers 6 and c are adapted for each spectrum to yield approximately equal S/N ratios for A,C, 
and D. This method allows calculation of a hardness ratio even when the source has an extremely steep or heavily absorbed 
spectrum. In practice, we use hardness ratios only when S/N> 4 in three bins {A, C, and D). Least-squares fits are used when 
the S/N ratio exceeds 5 in at least 6 bins. To determine the photon index F of the assumed power law from hardness ratios 
HRl and HR2, we fold through the PSPC response function a power law model with fixed Galactic absorption and derive the 
expected hardness ratios. We then vary the input T to minimize the sum of the quadratic differences between the measured 
and calculated hardness ratios, and quote errors corresponding to 68% confidence levels for a single parameter. A detailed 
description with tests of the method will be published in Schartel, Walter, & Fink (1995). More detailed tests can be found 
in Schartel (1994) 
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FIGURE CAPTIONS 



Figure 1. Photon count rate as a function of energy. As an example of our spectral stacking technique, we show the count 
rate for the total stacked spectrum (source plus background) of the 72 detected QSOs in the LBQS/RASS with 0.2 < « < 1. 
Vertical bars represent the errors in count rate, while horizontal bars show the energy range of the photons stacked in each 
energy bin. The area-normalized, stacked background spectrum for this sample is also shown. Poisson count rate errors for 
the stacked background are of similar size to the line thickness. 

Figure 2. Fits using 4 models to stacked X-ray spectra of 3 simulated QSO samples. EWA is best adapted to samples 
including non-detections, observed through different absorbing columns. The fractional difference spectra, (fit - input)/input, 
are shown in lower panels, a) Input spectra have a single power law slope (F = —2.5) but a range of Nh- EWA excels here, 
since the EWA fit is essentially identical to the input spectrum, b) Input spectra all have Nh = 2.5 x 10^° cm^^ but span 
a range of F. Fits using (Nh)^^ and (Nh)'^^ are identical to the input Nh in this case, and are omitted from the plot. 
The EWA model yields a superior fit to a model with free Nh. c) Input spectra have a range of both F and Nh. EWA is 
marginally superior to the model with counts- weighted (Nh)'^^ ■ However, the latter model can only be used when all sources 
are detected. 

Figure 3. X-ray photon index as a function of redshift. The best-fit (0.1 - 2.4 keV) mean power law photon index (F) 
from stacking is shown as a function of redshift for LBQS QSOs in 4 redshift ranges as observed in the RASS (filled circles 
with solid bars). The mean photon index for all 908 LBQS QSOs is also shown (filled circle with dashed bars). Vertical bars 
are the uncertainty in the photon index, while horizontal bars show the redshift range of the QSOs measured. 



